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We present a low-energy effective model for the charge degrees of freedom in two-leg t-J ladders. 
Starting from SU(2) mean-field theory, we exclude the spin degrees of freedom which have an energy 
gap. At low temperatures, the mean-field solution is the staggered-flux phase. For gapless charge 
excitations the effective theory is the Luther-Emery liquid. Our analysis is applicable at low doping 
and in the "physical" range of parameters t/J ~ 3 where there is only one massless mode in the 
charge sector and no massless modes in the spin sector. Within our model we make predictions about 
■yj " correlation exponents and the superconductivity order parameter, and discuss the comparison with 

■ the existing numerical results. 

< 
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Theoretical studies of t-J ladders have proven to be valuable for understanding Hi-T c cuprate compounds. While 
containing certain features of the two-dimensional t-J model, the ladders are quasi-one-dimensional, which greatly 
simplifies the treatment of the problem. Recently, cuprate compounds with ladder structures have been produced 
experimentally, and data on their electronic properties have been obtained Q| . The ladder compounds may be used for 
verification of models and mechanisms of superconductivity proposed for layered cuprates. The experimental works 
7-H , are now complemented by numerical results on t-J and Hubbard models PS. Together with experimental results, 
they provide a good testing ground for any analytical treatment. 

The challenge of theoretically solving the t-J model on ladders arises from its strongly correlated nature. In 
' the real ladder compounds the coupling is nearly isotropic, i.e. the interchain coupling parameters (t and J) are 
close to the intrachain coupling. Thus, the problem does not have a small parameter and cannot be treated by 
a standard perturbation theory. Several works exist based on starting from uncoupled chains and then including 
interchain hopping and spin exchange as perturbations While weak-coupling approach is the most consistent 

and controlled of the existing analytic methods, it is not completely reliable as the interchain coupling increases and 
approaches the single-chain bandwidth. We shall further comment on possible corrections to this treatment. 

In the present paper we employ the SU(2) slave-boson mean-field approach Hf?]]. Although not a controlled 
£j . approximation, we believe that it can correctly capture the low-energy physics of the systems with a spin gap. In 
the paper we specialize to two-leg ladders, but our treatment may be further extended to any even-leg ladders which 
are known to exhibit spin gap. The spin gap ensures that most of the fluctuations around the mean-field state are 
massive. The only massless fluctuation is the 1+1-dimcnsional abelian gauge field which can be explicitly included as 
a pair-binding potential. 

The general idea of slave-boson method is to represent the vacant sites (holes) by an auxiliary bosonic field, 
which allows us to rewrite the non-linear no-double-occupancy condition as a linear constraint in terms of fermions 
(representing spin degrees of freedom) and bosons (representing charge degrees of freedom) |^] . Introducing auxiliary 
bosons expands the Hilbert space of states, and the system aquires an additional gauge symmetry. The mean-field 
ansatz breaks this symmetry, which is restored for physical correlation functions after averaging over all gauge- 
equivalent configurations. 

We use the SU(2) version of the slave-boson construction developed earlier for the two-dimensional problem HQ. 
In this method, the auxiliary boson has two components (we call the corresponding degree of freedom isospin; it 
is distinct from the actual spin) which describe holes in the two different ways: either as sites with no fermions or 
as sites doubly occupied by fermions. Thus extended, the system has a SU{2) gauge symmetry (rotating isospin). 
We choose to use the SU(2) formalism instead of U(l) version of the slave-boson method developed previously for 
two-dimensional t-J model J§|| . We note that both analytic and numerical works H,||,|l^] point to a bipolaron picture 
where the holes are bound in pairs. As we shall see, this picture naturally emerges out of the 5/7(2) formulation as the 
confinement between two species of bosons, while the U(l) formalism fails to give the correct low-energy physics. We 
believe that the SU(2) mean-field theory has an advantage at low doping where it generalizes the SU{2) symmetry 
of t-J model at half-filling ||. A discussion of relation between U(l) and SU(2) approaches may be found in fij]. 
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We find that the low-temperature mean-field phase is the staggered-flux phase (similar to that found for the two- 
dimensional problem in 0J^]). The fluctuations about this mean-field state are described by gauge fields which have 
a gap, except for the in-phase fluctuations which form a 1+1-dimensional abelian gauge field. As a result, we find 
that the low-energy effective theory consists of two degenerate bands of holons with a short-range interaction and 
coupled to a U(l) gauge-field. The dynamics of this gauge field arises from its interactions with spinons. Spinons 
have a gap and, therefore, give a nonsingular dynamics to the gauge field with the energy scale J. The two bands of 
holons have opposite charges with respect to this gauge field and, therefore, form confined pairs. This leads us to the 
conclusion that the resulting theory for the hole excitations is the Luther-Emery liquid of hole pairs. This agrees with 
the bipolaronic picture of charge excitations proposed earlier both analytically and numerically |^,||,^0). It is known 
that the Luther-Emery liquid has two competing orders: superconducting singlet pairing (SS) and charge density 
wave (CDW) [p|J^Jll[|. We point out the necessity to distinguish between the hole density (two-particle operator) and 
the pair density (four-particle operator). While the product of the correlation exponents for the SS and CDW order 
parameters is equal to one when CDW is understood as hole-density correlations, this relation does not necessarily 
hold for four-particle pair-density correlations. The relation between the single- hole and pair CDW exponents depends 
on the degree of the overlap of the bipolaronic pairs. At low density of holes, when the pairs do not overlap, these two 
exponents coincide. On the other hand, in the limit of highly overlapping pairs we find that the effective exponents 
may differ by 2. This possibly explains the unexpected numerical results for the correlation exponents obtained by 
Noack et al. §]. 

Further, we discuss the possible implication of our model for the superconducting transition via pair condensation (of 
course, interladder correlations would be necessary). We describe the superconducting order by the nearest-neighbor 
order parameter Ay. We extend our discussion for a more general case of a weakly doped antiferromagnet on a 
bi-partite lattice with a spin gap. This class of systems includes all even-leg ladders as a particular case. We assume 
that at low temperature such a system is in the staggered flux phase, which results in two degenerate interacting 
holonic bands. Under these assupmtions we find that the order parameter obeys the modified d-wave relation: 



«A« = 0, (1) 

where the sum is performed over all nearest neighbors of a site i. This relation holds in the limit of zero doping, with 
corrections involving the hole concentration. This relation was first derived by S. C. Zhang as an exact result for the 
Hubbard model p2[ . Our derivation should be understood as a verification that the SU(2) mean-field approximation 
preserves this exact property. When specialized to the case of the two-leg ladder with isotropic coupling, the above 
equation becomes Aj_ = — 2Ay. This agrees very well with the earlier numerical results 

The rest of the paper is organized in three sections. In the first section, we review the SU (2) slave-boson method 
and present the results of the mean-field theory computations. In the second section we discuss the effective theory for 
the holons and the correlation exponents for SS and CDW pairing. Finally, the third part is devoted to the discussion 
of the modified d-wave relation for the superconducting order parameter. 

II. SU{2) MEAN-FIELD THEORY OF THE LADDER 

In this section we present the SU(2) mean-field theory for the t-J Hamiltonian 

H = J2 JiSiSj - - A n inj ) - tV(ci lCaj + h.c.)V (2) 

m 



on the ladder (Fig. |). 




FIG. 1. Two-leg t-J ladder. 



The sum is performed over nearest-neighbor site pairs {i,j}, c* ai and c Q ; are the electron creation and annihilation 
operators on site i (a is the spin index), Si is the electron spin, Si = \(y ai a a pCf3i, nt is the occupation nuber of the 
site i (m — c ai c a i), V is the projector onto the no-double-occupancy states (with rii < 1 for any i). t and J are the 
parameters of the Hamiltonian. In the real ladder compounds t and J are estimated to be about 4000K and 1300K 
respectively Jl]]. In what follows we assume that the interchain and the intrachain couplings are equal (t± = tii = t, 
J± = J|| = J) and as a realistic approximation we take t/ J = 3. 

Following the usual procedure of the SU(2) slave-boson method we introduce two fermionic and one bosonic 
isospin doublets on each site: 

*"-(£) **'(-%) (3) 

with the electronic operators written in terms of bosons hi and fermions ipai as 

1 



hU m (4) 



The resulting Hilbert space is larger than that of the original t-J system. To select the subspace of physical states 
(which is invariant under the Hamiltonian of the original system) we impose a linear constraint (replacing the non- 
linear no-double-occupancy constraint): 

(t^L^/w + h\Thi)\phys) = (5) 

(r are identical to the Pauli matrices er, but they act in the isospin space, and we denote them by a different letter to 
distinguish from the Pauli matrices a acting on true spin). On a given site, this constraint allows only three states: 
/i|0), /||0), and -^(b\ + b\flfl)\0) , which correspond to spin up, spin down electrons and a vacancy respectively. 
Thus in the SU(2) formulation, a vacancy may be represented by both two-spinon and no-spinon states corresponding 
to different isospins of the holon hi. 

Thus formulated, the extended system is invariant under an 577(2) gauge symmetry: 

Tpai i-> 9i1pai, hi h-> gihi. (6) 

This gauge symmetry acts on the isospin of fermions and bosons, and mixes creation operators f ai with the annihilation 
operators of opposite spin f- a i- 

Introducing the nearest-neighbor mean-field parameters 



U, 



= (^l)--Mh i h\) + i(h j ht) T l) (7) 



with 



1 

-1 



(8) 



and the Lagrange multipliers (fi = 1,2,3) for enforcing the linear constraint (||), the mean-field Hamiltonian 
becomes 



{ij} 



^Tr(C/ 4 ,4) + J -^ m U^ aj + ^hlUijhj + h.c.) 



+ E<(^» t ^« + / 4 t A)- (9) 



In the present paper we use combinatoric coefficients in (^) different from those used in |jM[l. Our choice of 
coefficients in (j^) gives the correct combinatoric factors for tadpole diagrams, and they differ from those obtained 
from the Hubbard-Stratonovich transformation. This ambiguity in the numeric factors of order one is of no practical 
significance, since it falls within the uncertainty of the mean-field approximation. 

The mean-field Hamiltonian (|^) is accurate only up to four-boson terms which we neglect in the usual way [H. 
The four-boson terms give rise to only short-range interactions between holons and can be omitted at this level of 
approximation. 

In Mv it was argued that the non-zero values of af correspond to Bose condensation. We do not expect Bose 
condensation in a quasi-one-dimensional system and set mean-field value af = 0. This implies that we in fact release 
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the no-double-occupancy constraint (it is satisfied only on average). Afterwards, the constraint may be imposed by 
including the fluctuations of the field af. 

Introducing the chemical potential fi, we arrive at the mean-field Hamiltonian 

H = J2 l-TriUijUl) + ^ m U^ aj + \{h\lJ l3 h 3 + h.c.)| - ^(h\hi ~ *)> ( 10 ) 
where 6 is the concentration of holes (doping). The matrices Uij have the form 

0y = (?V A ti)=u] i = ia ij G ij , (II) 



where ay are positive real numbers (amplitudes), Gij S SU(2) are 2x2 matrices. 

The mean-field saddle point at temperature T is found as the extremum of the free energy 

F[Un\ = -TlogTr^ exp(-H/T). (12) 

F[Uij] is invariant under SU(2) gauge transformations 

i * WiUijWj (13) 

for any set of SU(2) elements Wi. 

The mean-field solution breaks this gauge symmetry. Only gauge invariant quantitites correspond to physical 
observables. Any non-gauge-invariant expression will vanish after averaging over all gauge-equivalent configurations 
of Uij. 

Now we turn to describing possible phases. Phases should be parametrized by gauge-invariant functions of Uij. 
We assume that the translational symmetry is unbroken, i.e. a translation of the mean-field solution {Uij} along 
the ladder transforms it to a gauge equivalent configuration. We also assume that the symmetry of reflection about 
the ladder axis (interchanging the two legs) is also preserved in the mean-field solution. Under these assumptions 
all possible phases may be parametrized by four real parameters: the two amplitudes ay and a± (intrachain and 
interchain respectively), and two SU(2) order parameters: 



6 =l Tr ii G v> ( m ) 



(15) 



with the products taken along the contours Ti and T2 shown in Fig. [2] (the first product contains four matrices, the 
second one — eight matrices). 





-cM>- 

A 

FIG. 2. The two closed contours used in constructing mean-field order parameters. 

The meaning of the order parameter b is analogous to the cosine of the flux through plaquet in the U(l) formulation 
J. To explain this analogy we may introduce the 5/7(2) flux B defined by 




Ti 



(16) 



where the product starts and ends at a site A of the contour Ti (Fig. g). Then b = cos \ B\ (obviously, the direction 
of the vector B depends on the choice of the starting point A, but its magnitude \B\ does not). 
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The order parameter c measures the relative orientation of neighboring SU(2) fluxes. Namely, define the two fluxes 
B\ and B 2 through neighboring plaquets with a common starting point A (Fig. |J) . Then 

f / \ ( B B ) 

c= -Tr (exp^ •f)exp(-iB 2 • f)J = cos \B t \ cos \B 2 \ + } ' J sin \B X \ sin \B 2 \. (17) 
2 V / |i?i||i?2| 

Then we see that for any G^- 6 SU{2) the order parameters 6 and c are restricted to 

-1<&<1, (18) 



2b z 



1 < c < 1. 



(19) 



Thus, the SU(2) order may be represented by a point in a two-dimensional domain (Fig. ||) with the two corners 
representing the 7r-flux and the uRVB phases, the boundaries corresponding to the uniform flux (uF) and staggered 
flux (sF) phases analogous to their U(l) counterparts 0, but preserving translational and time-reversal symmetries. 
In the uF phase the neighboring SU (2) fluxes are parallel, in the sF phase they are antiparallel, and inside the shaded 
region in Fig. ^ they form angles ranging between and tt. 



uRVB 




FIG. 3. Space of mean-field phases. The two corners correspond to the 7r-flux and uRVB phases, boundaries — to the 
uniform-flux and staggered-flux phases. 

By numerically minimizing the free energy ( |l2| ) (at t/J = 3) we find the following mean- field phase diagram in the 
(6, T) coordinates: 



0.25J 




0.02 

FIG. 4. Mean-field phasee diagram at isotropic coupling and t/J — 3. 
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where 



1. denotes the high-temperature free spin phase (a = aj_ = 0); 

2. uRVB is the phase with an ^ 0, a± ^ 0, G L j = 1 (so that b = c= 1); 

3. D is the dimer phase: aj_ 7^ 0, an = 0; 

4. sF is the staggered-flux phase with a± ^ 0, an 7^ 0, — 1 < b < 1, c = 26 2 — 1. 

Numerically we find that the transition between D and sF phases is a very soft first-order transition. In fact, the 
dimer phase D has flat spectra for bosons and fermions and will be destroyed by the fluctuations (correlations along 
the ladder will appear). We shall disregard the dimer phase as an artifact of the mean-field approximation and for 
the rest of the paper we restrict our discussion to the sF phase. 

The sF phase may be described by different gauge-equivalent configurations of fTy. One of the translationally 
invariant configurations (analogous to d-wave pairing phase in the U(l) mean-field theory ||) is (Fig. ^): 

G± = in, G|| = i(cos |n + sin |r 2 ), (20) 

or, equivalcntly, 

u± = (« °A u, = ( » Hf). ( 21 ) 

\aj_ J 11 \a\\e '2 J 



G„ 

—0-^-0— 



Ql 1 ' 



—0-^-0— 

G„ 

FIG. 5. sF phase order parameters in the translationally invariant gauge. 
Further, we shall use a different, the so called "abelian" parametrization with 

G {1 =1, Gx = cos| + (-l) m *73sin| ) 



(Fig. ||) which corresponds to 

t/ll = 



Q|| 

an 



U± 



i(-i) m f 






-i(-i) m | 



(22) 



(23) 



where m is the number of the rung. This gauge fixing is not translationally invariant, but instead it has the property 
that all Uij commute. This choice of gauge resembles the staggered-flux U(l) phase ||. In the U(l) formalism, the 
staggered-flux and d-wave pairing phases are different, but their 5/7(2) counterparts are gauge-equivalent 0. 




FIG. 6. sF order parameters in the abelian gauge. 



6 



In the SU(2) mean-field theory bosonic and fermionic spectra are proportional to each other (with the scales t and 
J respectively). In the sF phase we obtain the spectrum 



El 
j 



El 
t 



±- 



(2ai|) 2 cos 2 fc + 2(2a|| 



)a± cos — cos k 



(24) 



(k G [0; 2ir] is the wave vector) and each of these two bands is doubly degenerate (four bands total). We remark that 
labeling the states by wave vectors k may depend on the gauge, but, because gauge-dependent shift of k involves 
equally bosons and fermions, the gauge-invariant quantities remain unchanged. 

The double degeneracy of states is a characteristic feature of the staggered-flux phase. It is due to the fact that 
a certain non-abelian subgroup of the full symmetry group of the extended Hamiltonian remains unbroken by the 
mean-field ansatz. In the staggered-flux phase the order parameters Uij are invariant under a U(l) subgroup of global 
isospin rotations (in the abelian gauge ( p2| ) - (23) these are simply global rotations by T3). Besides, there remains a 
particular symmetry of U^, which is the combination of time reversal (transforming fTy i— > U*j) and global isospin 
rotation (exchanging isospins up and down in the abelian gauge). This symmetry operation does not commute with 
the U(l) rotation, but extends them to a nonabelian group. The two degenerate bands form a two-dimensional 
representation of this group, with time reversal mapping one band onto the other. 

The typical numerical values for ay, a± and cos ^ are ay = 0.5, a± = 0.9, cos ^ = 0.4 (fonud by minimizing free 
energy at 8 = 0.05, T = 0.1J, t/J = 3). This means that the upper bands of the spectrum are separated from the 
lower bands by a gap of order J for fermions and of order t for bosons (Fig. [7J). The fermionic spectrum is half- filled, 
i.e. the lower bands are completely filled whereas the upper bands are empty. The fermionic excitations have a gap of 
order of J. This agrees with the prediction of spin gap in two-leg ladders and is crucial for the stability of the phase. 




FIG. 7. Typical fermionic-bosonic spectrum. Both bands are doubly degenerate. 



Finally, we need to include the fluctuations about the mean-field phase. These fluctuations are described by the 
spatial SU(2) gauge fields Ay = U^ SUij defined on the links and their temporal counterparts a, = a£V M defined 
on the sites H The temporal components of the gauge field a, coincide with the constraint-fixing Lagrange 
multipliers in (|9|) when expressed in units of J. The effective action for the fields Ay and a, is invariant under the 
(time-dependent) gauge transformations 

Aij 1 ^ Wii^AijWjit), a, >-> a, + id t W t {t) (25) 

for arbitrary SU(2) matrices Wi(t) defined on sites. 

In the sF phase the SU(2) gauge symmetry is broken down to U(l) global symmetry. In the abelian gauge, this 
residual symmetry is realized by rotations by T3. As shown in the Appendix B of |Q, the free energy will contain 
terms proportional to Tr Y[ Uij , where the products are taken along closed loops on the lattice. When expanded in 
gauge-field fluctuations, these terms give rise to mass for the gauge modes proportional to t\ and T2- Thus, massless 
modes of the gauge field must be proportional to T3 M, and we may treat Aj and a% as U(l) gauge fields. 
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To proceed further, we may make use of the symmetry (25) (now an £7(1) gauge symmetry with all Wi(t) being 
rotations by T3). The analysis appears to be particularly simple in the gauge where Aij — across the rungs (Fig. ||). 



^ A ( ° & 
1 y ^ J 



a 



(n) A (n) 



FIG. 8. Gauge fields describing the fluctuations of the order parameter. We choose the gauge with A± — 0. 



Then the remaining gauge fields split into the in-phase and out-of-phase modes 

A ± =4>±4>\ a ± = ^±^ (26) 

(where the superscripts (I) and (77) label the two legs of the ladder). The dynamics of the gauge fields arises from 
their coupling to the spinons and can be found by computing the polarization diagram (Fig. ||), solid lines denote 
spinons): 




FIG. 9. Lowest-order diagram responsible for the dynamics of the gauge field. Solid lines denote spinon propagators. This 
diagram gives mass to the out-of phase mode and 1+1-dimensional QED dynamics for the in-phase mode. 

In contrast to the two-dimensional model, where there exist massless transverse fluctuations of the gauge field, it is 
not the case in the ladder. The ladder geometry restricts the transverse wavevector to two values k± = and = 7r 
which correspond to the modes (A + ,a + ) and (A_,a_) respectively. The out-of-phase modes A- and a_ (describing 
the fluctuations of the flux ip through the plaquet) acquire a finite mass of order J and, therefore, can be neglected, 
giving only a short-range interactions between holons. 

On the other hand, the modes A + and a+ become 1+1 QED gauge fields with the long- wavelength action 

S = T* j + dxa+)2 dx dt (27) 

with J* of order J. This action arises from expanding the polarization diagram in Fig. [)J The particular value of J* 
depends on the values of mean-field parameters, and for the typical values cited above differs from J only by a factor 
of order unity. 

Since the two lower bosonic bands have different isospin, they have opposite charges with respect to the gauge field 
(A-|_,a-|_). In 1+1 dimension, electromagnetic field leads to a confining (linearly growing with the distance) potential 
between charges: 

U(r) = ±J*\r\. (28) 

In the limit of low hole density, the bosons will form isospin-neutral dipole pairs (bipolarons) which produce no field 
outside each pair. Therefore, bipolarons will interact only by short-range forces. 

We may estimate the size of bipolaron by solving a simple quantum-mechanical problem of two particles interacting 
via the potential (|2|). The holon hopping amplitude is of order i, therefore the kinetic energy of bipolaron is of 
order t^~^ ir . The potential energy, on the other hand, is of order J£, pa ir- Thus from variational principle we find 
that the size of bipolaron is S, pa i r ~ (t/J) 1 ' 3 up to a factor of order one. For our assumption t/J = 3, this gives 
£,pair ~ 1, which means that now the other short-range (repulsive) terms, which we omitted before, give a comparable 
contribution. Due to the no-double-occupancy constraint, the two species of bosons must be subject to a substantial 
on-site repulsion. It increases the estimated size of the pair by several lattice spacings, but does not change the 
long-range attractive force. Since we omitted the short-range part of the interaction from the very beginning, we 
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cannot compute the size of the pair more precisely. From the numerics g we know that the characteristic decay 
length for the pairing correlation function is about four lattice spacings, in agreement with our discussion. 

A simple classical explanation of the confinement may be obtained from the picture of fluctuating singlet bonds 
(somewhat in the spirit of JTo|). We may think of the spin structure of the Heisenberg antiferromagnet on the ladder 
as of spins forming fluctuating singlet bonds with nearest neighbors (Fig. [To| ). 




FIG. 10. Fluctuating singlets in the spin ladder. 



Once we have a single hole, it leads to the appearance of a localized spin, which costs a finite energy of order 
J. When putting two holes, the spins between the holes must form singlet bonds in a non-favorable way without a 
freedom to fluctuate (Fig. |ll|). 




FIG. 11. Misplaced singlet bonds between two holes in the t-J ladder. The energy cost is proportional to the distance 
between the holes. 

Naturally, this string of singlets costs certain energy J*r, where r is the distance between the holes and J* is of 
order J. We believe that this naive picture gives a correct understanding of the holon confinement which we derived 
starting from the mean-field sF phase. 



III. LUTHER-EMERY LIQUID AND CORRELATION EXPONENTS 

In the previous section we have shown that the low-energy excitations in our model are pairs of holons — bipolarons 
— which are bound by long-range confining interaction. Due to complete screening of the confining interaction by 
a single particle in one dimension, bipolarons interact only via short-range forces, and this recovers the picture of 
Luther-Emery liquid of hole pairs proposed earlier in |fl0|| . 

In fact, this simple picture is valid at low hole concentration S <C £,pair> when bipolarons do not overlap. The 
purpose of this section is to show that as overlap increases, the charge-density- wave exponents for single-particle and 
for pair densities may differ. 

For quantitative description of bipolarons we introduce the three different correlation exponents as follows. Let 
n(x) = + n^ 77 ) — 26 be the fluctuation of the number of holes on the rung (x is the coordinate along the ladder, 
superscripts refer to the two chains). Let further n pa ir — {nW — (5)(n^ // ' — 8) be the fluctuation of the probability 
that both sites of the same rung contain holes. Finally, let A(x) = c'pc^ 1 " 1 — c^'c^ 1 ' be the singlet superconducting 
order parameter on the rung at the position x. Then define the correlation exponents a\, a 2 and 7 by 

(n(x)n(y)) 2kF *\x~y\- a \ (29) 
impair {p^jf^pair (y))2k F oc \x-y\ " 2 , (30) 

(ti(x)A(y))cx\x-y\-\ (31) 

where (. . -)2k F is the coefficient at coa[2kp(x — y)] in the expansion of the correlation function Q]: 

(n{x)n(y)) = A \x - + A x \x - y\~ 8 ^ cos[2k F {x - y)] + 

+A 2 \x - y\- 02 cos[4fc F (a; -y)} + ... (32) 
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The correlation function (|29|) describes charge-density-wave (CDW) ordering, and the correlation function ( J3l| ) - 
singlet superconductivity (SS) ordering. 

By kp we denote the "Fermi wavevector" for the holons in the lower bands, so that 2kp = 2ir5. For comparison 
with other works, we must remark that 2kp in our notation corresponds to 4fcp in the notation of S and S. 

In the limit of dilute gas of bipolarons (the doping is much less than the inverse size of a bipolaron: S <C £p a i r ) we 
may describe the low-energy states of the system in terms of bipolaron creation and annihilation operators (similarly 
to the large \U\ limit in the attractive Hubbard model |ll],[l3)). Bipolarons interact repulsively and, from a naive 
classical picture, it is likely that the repulsion is nearly hard-core. In this limit a\ — a 2 , because on the bipolaronic 
subspace of the total Hilbert space the matrix elements of the operators n{x) and n pa i r (x) differ only by a numerical 
factor (the probability of the two holons in the pair to occupy the same rung). For hard-core repulsion = 2, since 
hard-core bosons can be mapped to free fermions by a Jordan- Wigner-type transformation, so that density-density 
correlations coincide with those of free fermions. From the theory of Luther-Emery liquid it is known that CDW 
and SS correlations are described by dual phases, and the corresponding exponents are therefore reciprocal: 7 = 1/a 
jll],[l3| . A more detailed discussion of the Luther- Emery theory of bipolaronic excitations may be found in jlo| . 

The opposite limit of highly overlapping bipolarons is more subtle. In this limit (S <C (i/J) 1 / 3 ) bipolarons can 
exchange particles. Exchanging a particle would cost interaction energy of order Ji5 _1 , while the gain of kinetic energy 
would be of order tS 2 . Thus, in this limit we cannot speak of isolated bipolarons, but rather of two species of bosons 
with attraction much smaller than the bandwidth. We suggest that in this case we may replace the screened long-range 
interaction by a short-range one. In the long-wavelength limit the system may be described as two Luther-Emery 
liquids with a weak attraction: 

H = Hi + H 2 + Hi n t , (33) 
Hi = i- j dxlvjiVifii) 2 +v N {Ve t f], i = 1,2, (34) 



Hint 



-V J dxpi(x)p 2 (x), (35) 



where ipi and 9i are dual phases, 



7T 

mi x )> 9j(y)] = S l:j i-sign(x - y), (36) 



vj and vn are the parameters depending on the short-range properties of the interaction, and pi (x) are the density 
fluctuations expressed by Q 

Pi(x) — —\76i(x) + 2p cos{2kpx + 29i(x)) + higher order terms. (37) 

The interaction ( |35| ) contains a term proportional to cos 2 (#1 (a;) — #2(2:)) and locks the relative phase 9- = 9\ — 9 2 - 
The only remaining gapless mode is the in-phase fluctuations {9 + = 9\ + 9 2 and ip + — ipi + ip 2 ) corresponding ot 
propagation of bipolarons. 

If we at first approximation neglect higher-order terms in the density expansion ( |37| ) , we find 

(AHx)A(y)) - ( e %+ (*)-?+(!/)]}, ( 38 ) 
(n(x)n(y)) 2kF ~ (e^+W-Mv)]), (39) 

(n pai r(x)n pair (y)) 2 k F ~ (\79 + (x)^9 + (y)e^^~ e +^), (40) 

which yields a 2 — ct\ + 2, 7 = l/ot\. This would explain the numerical results of Q who found a 2 ~ 7 ~ 2 (at 
S = 1/8, ^,p a i r ~ 4). Our prediction of power-law correlations (29) also explains the small, but relatively narrow peak 
at 2kp in the Fourier transform of (n(x)n(y)) in [H. 

As it was pointed out by Haldane H] , in general one must also include higher-order terms in the density expansion 
(|37|). These terms proportional to cos[m(2kpx + 29i(x))] with m > 1 are absent in the free fermion theory (and, 
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consequently, in hard-core boson theory), but arise as we include interactions mixing left- and right-moving excitations. 
They produce terms proportional to exp i[m9i(x) — (m± 1)02 (a:)] in the n pa i r (x)2k F expansion, and give a contribution 
to (n pa i r (x)n pa i r (y)}2k F decaying with the exponent o.\ instead of a\ + 2. This effect that originally higher-order 
terms result in leading correlation exponents is not paradoxical in view of the crossover to the dilute limit where the 
binding interaction is strong and the higher-order corrections to ( |37| ) play a dominating role, we suggest that the 
crossover from ai = o.\ in the dilute limit to 0,2 = ol\ + 2 in the weak-coupling limitis governed by the overlap of 
bipolarons. Namely, 

A B 

(n pair (x)n pair (y))2 kF = {x _ y)ai + (x _^ )ai+2 ( 41 ) 

with relative weights of A and B depending on the average pair overlap (A ^> B at £, pa ir3 "C 1 and A <C B at 
£,paird 3* 1). The actual behavior of the coefficients A and B strongly depends on the short-scale features of the 
interaction, and cannot be found in our rude treatment. 

The whole discussion of this section is equally applicable to the negative-U (attractive) Hubbard model. In the 
low-density (or large \U\) limit the exponents ol\ and U2 coincide, while in the small U limit we expect a crossover 
( ftl| ) to «2 = Oil + 2. In other words, because of screening, the long-range gauge interaction between holons in one 
dimension leads to the same behavior at large distances as a short-range attraction. 

Finally, we comment on our disagreement with the prediction of Nagaosa || that the correlations 
(n P air(x)n pa ir(y)}2k F decay exponentially (note again that 2kp in our notation corresponds to Akp in the notation of 
||). The disagreement my be explained from the fact that Nagaosa starts from two uncoupled chains and treats the 
interchain couplings (t± and Jj_) as perturbations. In that picture, power-law correlations (n pa i r (x)n pa ir{y)) 2k F will 
appear as a correction for the nonlinearity of the spectrum as the coupling increases and approaches the bandwidth. In 
contrast to the weak-coupling approach, our model starts directly from diagonalizing a strong-coupling Hamiltonian 
(|lo|), and the correlations (^) are present from the very beginning. 



IV. MODIFIED D-WAVE RELATION ON SUPERCONDUCTING ORDER PARAMETER 

In this section we verify that our approximation scheme is consistent with the exact relation for the superconducting 
order parameter derived by S. C. Zhang for the Hubbard model jlj] (and later translated to t-J model in 
Let us define the pairing operator 

Ay = a n c n - CifCjj. (42) 

and consider the quantity 

Ag> = <2|Ay|0>, (43) 

where |0) and |2) are the ground states with zero and two holes respectively. Should a superconducting transition 
happen, A^y will become the superconducting order parameter. 

From now on, we restrict i and j to be nearest-neighbor sites. Zhang's results states that on a bipartite lattice in 
the limit of zero doping obey the relation 

E4 0) = °> ( 44 ) 

3 

where the sum is over the nearest neighbors of the site i. On the two-dimensional square lattice this implies the 
d-wave symmetry of pairing; thus we may call Eq.(f44|) the modified d-wave relation. 

Below we rederive this result within the SU(2) slave-boson mean-field approximation. For the sake of generality, 
we extend our further discussion to the t-J model on an arbitrary quasi-one-dimensional bi-partite lattice, provided 
it exhibits a spin gap (the most popular examples of this type are even- leg ladders). Further, assume that the low- 
temperature mean-field phase is analogous to the sF phase of the two-leg ladder. Namely, we require that the SU (2) 
order parameter may be brought to the diagonal form (by a suitable choice of gauge): 

Uij=( X ^ (45) 



-x 



i j 
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This requirement means that spinons and holons form doubly degenerate bands related by the symmetry of simul- 
taneous time-reversal and isospin flip. In such a phase a superconductivity may evolve by pair formation between 
the holons at the bottom of the two lowest bands (of course, superconductivity is possible only when stabilized by 
inter-ladder interactions, see e.g. 

For simplicity, let the coupling be isotropic (t and J are the same on all links), as it was assumed in the previous 
sections. At the end of this section we shall extend the result to non-isotropic coupling. 

Using our slave-boson representation, we express A^- (for nearest-neigbor i and j) in terms of spinons and holons 

as 



\ [{h\^ 2i ){h)^) - (/4Vii)(^)] • (46) 
At low doping, the fermionic part of the correlation function (4.2) may be replaced by the mean-field order parameters 

Xij = (hiflj + hiflj), (47) 

and we find 



A; 



\ (xiMA&ft) + X^2|4A-|0>) > (48) 



where |0) and |2) now denote the states in the holonic sector. Let b\(k) and &2(fc) be the operators destroying holons 
at a wave vector k in the two lowest bands (subscripts denote the isospin). Then the single-pair wave- function |2) has 
the form 

|2> = J ^o(k)b{(k)bl(-k)\0), (49) 
where ^o(k) is the relative wave function of the two holons in a pair. Eq. (f48|) becomes 

Ag } = J ^*o(A)Ay(fc), (50) 

where 

Ay(fc) - \ (x«(0|6 2 (-fc)6i(fc)6 t 1 4j|0)+X^(0|6 2 (-*)6i(fc)46 t M |0)) = 
= \ (x«(0|62(-fc)6j i |0)(0|6 1 (fc)6l i |0)+ X * i (0|6 2 (-fc)4|0)(0|6i(fc)6l J -|0)) . (51) 

Since bi(k) and &2(fc) are related by the time-reversal symmetry (accompanied by a gauge transformation), 
(0| 62 ( — fc)^2i|0) = (~ l) 4 (0|^ii^i(^)|0), and for nearest-neghbor sites i and j 

Atf(fc) - \ ((-iyx«(0|&y& t 1 (fc)|0)(0|6 1 (fc)& t li |0) + (-l) i X * i (0|6 M 6 t 1 (fc)|0)(0|6 1 (fc)6 t lj .|0)) 

= (-m^i-xiAihi+xtAMW' ( 52 ) 

where \k) is a single- holon plane wave created by b\(k). The state \k) is an eigenvector of the free Hamiltonian 
proportional to the bosonic part of Eq. ( jToj ) 

H« = \ E( 6 ky& U + bljxljhi) (53) 

(with the sum performed over nearest-neighbor pairs of sites) . Therefore 

^A^(fc) = (-iy(k\[blb u ,H ]\k) =0, (54) 
j 

where the sum is over the nearest neighbors of the site i. This immediately implies the result (H). 
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The above derivation holds also at a non-zero temperature (with ground-state expectation values replaced by 



thermal averages). However, it is strictly limited to zero doping: at finite doping Eq.(47) is no longer valid. A 
remarkable feature of the relation (Q) is its independence of the coupling parameters t and J. 
For the two-leg ladder the modified d-wave relation (Q) turns into 

A_L = -2Ay. (55) 

This result agrees with the available numerical results J2|,[5| . It would provide a good test for possible numerical models 
on ladders with a higher number of legs. 

The d-wave relation ( fl4| ) may be easily extended to a non-isotropic coupling. In fact, the whole slave boson mean- 
field theory may be rederived for arbitrary coupling constants t%j and Jy , differing at different links. One just needs 
to replace t and J in Eqs. (^) - ([To]) by Uj and Jy. The whole argument of this section may be repeated downto 
Eq.(p3[) which we must now replace by the Hamiltonian 

H « = \Y. *« i^iXijhj + b^jbu). (56) 
Finally, this leads to the following generalization of the d-wave relation (Q) : 

E*X 0) = - ( 57 ) 

3 

This equation implies that as the hopping on a link increases, the weight of the superconducting oreder parameter 
on this link decreases. Of course, Eq.(|57]) may also be derived exactly for the Hubbard model (with the site- and 
link-dependent U and t) by the method of [Ej|. 



V. CONCLUSION 



We presented the low-energy effective theory for charge excitations in two-leg t-J ladder, based on the mean-field 
treatment of spin degrees of freedom. We found that the SU (2) slave-boson formalism predicts bipolaronic picture 
of charge excitations, as expected from earlier analytic and numerical works. While capturing well the low-energy 
physics, our approximation is not reliable for spin and single-electron excitations which have a gap of order J. In 
the framework of the model developed in the paper, single-hole excitations may be constructed as holon-spinon pairs 
bound by confining gauge-field interaction. The similarity of the boson-fermion spectrum in Fig. ^ to the single-hole 
spectrum found numerically in |l(J makes this possibility very appealing. However, spin structure of the ladder was 
included in the lowest mean-field order, and to restore correctly spin excitations requires a more elaborate treatment. 
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